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impacts  of  remote  forcing,  model  resolution  and  bathymetry  on 
current  predictions  at  two  moorings  located  on  the  shelf  of  the 
Monterey  Bay  area  are  Investigated.  We  consider  three  Mon¬ 
terey  Bay  model  configurations  which  differ  in  resolution  and 
bathymetry  representation,  and  we  specify  open  boundary  condi¬ 
tions  for  these  three  configurations  from  two  larger  scale  models, 
which  have  different  accuracy  in  the  representation  of  the  remote 
forcing  (in  the  form  of  poleward  propagating  along  the  coast 
coastally-trapped  Kelvin  type  waves). 

Comparisons  of  correlations  between  observed  and  model 
currents  as  well  as  visual  comparisons  show  that  the  most  critical 
element  in  reproducing  currents  on  the  shelf  is  accurate  represen¬ 
tation  of  the  remote  forcing.  Our  results  also  show  that  accurate 
representation  of  bathymetry  is  the  second  most  critical  factor  in 
reproducing  observed  currents. 

Published  by  Elsevier  B.V. 


1.  Introduction 

The  objectives  of  the  August  2006  field  experiment,  called  Adaptive  Sampling  and  Prediction 
(ASAP),  were  mostly  focused  on  the  study  of  the  properties  of  the  upwelling  center  at  Afio  Nuevo 
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Fig.  1.  (A)  The  Monterey  Bay  modeling  domain  with  locations  of  ADCP1,  ADCP2  and  the  ASAP  glider  sampling  domain  to  the 
north  of  the  ADCPs.  (B)  Grid  resolution  around  ADCPs  for  MBS1  and  MBS2  configurations.  (Q  Grid  resolution  around  ADCPs  for 
MBS3. 

to  the  north  of  the  Monterey  Bay  (Ramp  et  aL,  2011;  Leonard  et  al„  2010;  Shulman  et  al.t  2010).  For 
this  reason,  the  extensive  sampling  was  conducted  inside  of  an  approximately  1000  km2  box  (Fig.  1), 
where  a  fleet  of  ten  gliders  under  autonomous  control  were  deployed  for  a  period  of  30  days,  and 
research  aircraft  observed  the  fluxes  through  the  sea  surface.  Two  bottom-mounted  acoustic  Doppler 
current  profilers  (ADCPs)  were  deployed  about  6.5  km  apart  to  the  south  of  the  ASAP  box  (Fig.  1 ),  to 
monitor  the  currents  over  the  continental  shelf  (Ramp  et  al„  2011).  Despite  being  only  6.5  km  apart, 
ADCP  moorings  1  and  2  responded  differently  to  the  sequence  of  upwelling  favorable  winds  separated 
by  brief  relaxations.  Predictions  from  three  simulations  of  the  Monterey  Bay  area  based  on  the  Har¬ 
vard  Ocean  Prediction  System  (HOPS),  the  Regional  Ocean  Modeling  System  (ROMS),  and  the  Navy 
Coastal  Ocean  Model  (NCOM)  were  quantitatively  compared  with  the  observed  currents  at  the  two 
moorings’  locations  on  the  shelf  (Ramp  et  al„  201 1 ).  All  three  model  simulations,  with  well-established 
performance  at  larger  space  and  time  scales,  had  difficulty  reproducing  the  current  variability  in  this 
small  sample  region  around  mooring  locations  (with  relatively  better  performance  in  the  alongshore 
than  a  cross-shore  directions).  It  was  speculated  that  one  of  the  reasons  is  that  model  open  bound¬ 
ary  conditions  could  not  capture  remote  forcing  in  the  form  of  alongshore  pressure  gradient  forces 
or  coastally-trapped  waves,  which  propagate  from  south  to  north  with  the  coast  on  the  right  in  this 
region.  Other  considered  reasons  were  that  very  high  horizontal  resolution  (at  least  0.5  km)  and  more 
accurate  representation  of  bathymetry  are  needed  to  reproduce  currents  variability  on  the  continen¬ 
tal  shelf  (based  on  relatively  better  performance  of  the  finer  resolution  simulation  based  on  the  HOPS 
system).  Because  the  three  considered  modeling  systems  had  so  many  differences  in  specification  of 
open  boundary  conditions,  data  assimilation  schemes,  bathymetry,  horizontal  and  vertical  resolution, 
parameterization,  and  ways  of  applying  atmospheric  conditions,  it  was  difficult  to  sort  out  reasons  for 
model  difficulties  in  reproducing  currents  in  the  framework  of  the  Ramp  et  at.  (201 1 )  study. 

The  objective  of  this  short  follow  up  to  Ramp  et  al.  (2011)  paper  is  to  use  one  model  (NCOM) 
and  identify  the  impacts  of  remote  forcing,  model  resolution  and  bathymetry  representation  on  the 
model  current  predictions  on  the  shelf.  For  doing  these  we  consider  three  NCOM  configurations  which 
differ  in  resolution  and  bathymetry  representation  and  we  specify  open  boundary  conditions  for  these 
three  configurations  from  two  different  larger  scale  models,  which  have  different  accuracy  in  the 
representation  of  the  remote  forcing  (the  coastally-trapped  Kelvin  type  waves). 

The  structure  of  the  paper  is  as  follows:  Section  2  describes  observations,  Monterey  Bay  model 
configurations  and  open  boundary  conditions  used  in  this  study.  The  design  of  model  runs  is  described 
in  Section  3,  Section  4  presents  results  of  experiments  and  Section  5  is  devoted  to  conclusions  and 
discussions. 


I.  Shulman  etal. /  Dynamics  of  Atmospheres  and  Oceans  6 1 -62  (2013)  35-45 


37 


m 


I22.5W  122W  122.5W  122W  I22.5W  122W 

Fig.  2.  Bathymetry  for  model  configurations  MBS1 -MBS3. 


2.  Methods 

2.1.  ADCP  data 

A  comprehensive  description  of  ADCP  data  is  presented  in  Ramp  et  al.  (2011).  Here  we  provide 
a  brief  of  this  description.  The  moorings  deployed  for  the  August  2006  experiment  consisted  of  two 
300  kHz  acoustic  Doppler  current  profilers  (ADCPs)  in  trawl-resistant  bottom  mounts  (TRBMs)  located 
6.55  km  apart.  ADCP  1  was  located  at  36°  55.336' N,  122°  07.344' W.  54  m  depth,  and  ADCP  2  was 
located  at  36°  53.683'  N,  122°  11.244  W,  92  m  depth  (Fig.  1).  The  data  were  sampled  in  4-m  bins  from 
the  bottom  to  the  surface;  however  the  uppermost  3  bins  (12  m)  were  lost  due  to  side-lobe  reflection. 
The  basic  1-min  time  series  were  obtained  by  ensemble  averaging  60  1-s  pings  in  the  instrument. 
Subsequent  averaging  and  filtering  to  separate  the  frequency  bands  of  interest  were  accomplished 
in  the  laboratory.  An  inverse-Fourier  truncation  method  was  used  to  remove  internal  waves  with 
periods  of  five  hours  or  less.  The  residual  time  series  were  filtered  again  with  the  half-power  point  at 
33  h  (Beardsley  etal.,  1985)  to  separate  the  tidal  and  inertial  motions  from  the  wind  forced,  mesoscale, 
and  mean  currents. 

22.  The  NCOM  Monterey  Bay  model  configurations 

The  NCOM  is  a  primitive-equation,  3D,  hydrostatic  model.  It  uses  the  Mellor-Yamada  level  2.5 
turbulence  closure  scheme,  and  the  Smagorinsky  formulation  for  horizontal  mixing  (Martin,  2000). 

Three  NCOM  configurations  of  the  Monterey  Bay  area  are  used  in  this  study.  All  three  set  ups  cover 
the  same  domain  (shown  in  Fig.  1 ).  They  are  all  forced  with  surface  fluxes  from  the  Coupled  Ocean  and 
Atmospheric  Mesoscale  Prediction  System  (COAMPS)  (Doyle  et  al.,  2009)  atmospheric  model  at  3  km 
horizontal  resolution.  The  bathymetry  with  the  horizontal  resolution  of  approximately  0.01°  (~1  km) 
was  interpolated  to  the  centers  of  the  horizontal  grids  to  derive  depths  for  all  of  the  described  below 
three  model  configurations. 

2.2. 1.  First  Monterey  Bay  model  configuration  (MBS1) 

The  first  configuration  is  the  model  set  up  used  in  Ramp  et  al.  (201 1 ).  This  configuration  is  noted 
as  MBS1.  The  MBS1  has  a  curvilinear  orthogonal  horizontal  grid  with  resolution  ranging  from  1  to 
4  km.  near  the  ADPCs  the  resolution  is  about  1.5  km  x  2.5  km  (Fig.  IB).  There  are  30  sigma-coordinate 
vertical  levels.  Because  of  the  sigma-coordinate  vertical  system,  the  bathymetry  has  to  be  smoothed 
to  minimize  problems  associated  with  sigma-coordinates  over  steep  topography  (Haney,  1991).  The 
bathymetry  for  the  MBS1  is  shown  in  Fig.  2. 

222.  Second  Monterey  Bay  model  configuration  (MBS2) 

The  second  configuration  (noted  by  MBS2)  has  the  same  horizontal  resolution  as  the  MBS1 .  but  with 
a  different  vertical  coordinate  system.  MBS2  has  40  hybrid  (sigma-z)  levels  total.  1 9  sigma  layers  from 
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the  surface  down  to  138-m  depth  and  21  z-levels  between  138  m  and  bottom.  In  this  case,  with  z  levels 
on  the  bottom,  there  is  no  need  to  smooth  bathymetry,  and  bathymetry  used  for  MBS2  configuration 
is  shown  in  Fig.  2. 

22.3.  Third  Monterey  Bay  model  configuration  (MBS3) 

The  third  configuration  (noted  by  MBS3)  has  a  finer  resolution  horizontal  grid  in  comparison  to 
MBS1  and  MBS2.  The  grid  is  also  a  curvilinear  orthogonal  grid  with  a  resolution  of  0.4  km  x  1.5  km, 
and  near  the  ADPCs  the  resolution  is  about  0.5  km  x  0.7  km  (Fig.  1C).  The  vertical  grid  consists  of  49 
hybrid  (sigma-z)  levels  total,  24  sigma  layers  from  the  surface  down  to  146-m  depth  and  25  z-levels 
between  146  m  and  bottom.  The  bathymetry  for  MBS3  is  shown  in  Fig.  2. 

2.3.  Open  boundary  data  and  remote  forcing 

Two  larger  scale  model’s  outputs  are  used  to  provide  boundary  conditions  for  the  Monterey  Bay 
model  configurations  described  in  Section  2.2.  The  first  model  is  the  NCOM-based  regional  model  of 
the  California  Current  (NCOM  CCS,  Shulman  et  al.,  2007).  The  NCOM  CCS  has  a  horizontal  resolution 
of  about  9  km  and  uses  the  same  40-layer,  hybrid  vertical  grid  used  in  the  MBS2  set  up.  The  NCOM 
CCS  domain  extends  from  30  N  to  49°  N  of  latitude  and  from  the  coast  to  1 35°  W  of  longitude.  The 
NCOM  CCS  model  is  forced  with  atmospheric  products  derived  from  the  COAMPS  predictions  (Doyle 
et  al„  2009).  The  NCOM  CCS  assimilates  three-dimensional  temperature  and  salinity  observations 
derived  from  the  Modular  Ocean  Data  Assimilation  System  (MODAS;  Fox  et  al.,  2002).  Open  boundary 
conditions  for  the  NCOM  CCS  model  are  derived  from  global  NCOM  model  (Rhodes  etaL2002;  Barron 
et  al.,  2004). 

The  second  larger-scale  model  is  the  global  model  based  on  the  Hybrid  Coordinate  Ocean  Model 
(HYCOM)  (Chassignet  et  al.,  2009;  Metzger  et  al.,  2010).  The  model  has  ~1/12°  horizontal  resolution 
and  32  vertical  layers.  It  uses  the  Navy  Coupled  Ocean  Data  Assimilation  (NCODA)  system  for  assim¬ 
ilation  of  satellite  surface  temperature,  as  well  as  available  in  situ  vertical  temperature  and  salinity 
profiles  from  XBTs,  Argo  floats,  moored  buoys  and  gliders  from  the  Global  Ocean  Data  Assimilation 
Experiment  (GODAE)  data  set  (Cummings,  2005;  Cummings  et  al.,  2009).  The  model  is  forced  with 
atmospheric  fluxes  from  the  Navy  Global  Atmospheric  Prediction  System  (NOGAPS;  Rosmond  et  al., 
2002). 

In  Ramp  et  al.  (2011,  Fig.  12),  the  comparison  of  sea  surface  height  (SSH)  anomalies  along  the 
California  coast  estimated  from  the  NCOM  CCS  and  the  global  HYCOM  with  the  coastal  SSH  observa¬ 
tions  for  May-September  of  2006  are  presented.  Both  models  and  observations  show  strong  positive 
SSH  anomalies  during  August  2006  (Ramp  et  al.,  201 1 ;  Shulman  et  al.,  2010).  Both  models  also  show 
that  these  anomalies  were  propagating  poleward  at  1-3 ms-1,  which  is  in  agreement  with  previous 
estimates  of  phase  seed  for  the  first  baroclinic  coastally-trapped  Kelvin  wave  (e.g.  Chelton  and  Davis, 
1982;  Denbo  and  Allen,  1987;  Spillane  et  al.,  1987;  Zamudio  et  al.,  2008,  2011).  These  anomalies 
therefore  represent  remote  forcing  which  will  impact  the  Monterey  Bay  area,  and  a  smaller  region 
around  ADCP 1  and  2  on  the  shelf.  Comparing  the  model  SSH  with  the  observed  SSH  at  Monterey  and 
San  Diego  tide  stations  shows  that  the  NCOM  CCS  underestimated  these  propagating  anomalies  with 
respect  to  both  the  observations  and  HYCOM,  which  closely-tracked  the  observed  fluctuations  (see 
Fig.  3  which  is  Fig.  1 2  of  Ramp  et  al..  201 1 ,  and  reproduced  here  for  clarity  presentation).  Therefore,  the 
use  of  the  HYCOM  output  as  open  boundary  conditions  for  the  Monterey  Bay  model  (configurations 
MBS1-3  in  Section  2.2)  will  provide  more  accurate  representation  of  remote  forcing  on  the  southern 
open  boundary  in  comparison  to  using  the  NCOM  CCS. 

Sea  surface  elevations  and  vertically  averaged  velocities  (barotropic  information)  from  the  Mon¬ 
terey  Bay  model  and  larger  scale  model  (NCOM  CCS  or  HYCOM)  are  coupled  through  the  Flather(1976) 
boundary  condition: 


(1) 


where  r;  and  un  are  the  Monterey  Bay  model  (configurations  MSI -3,  Section  2.2)  sea  surface  height 
and  the  vertically  averaged  normal  component  of  velocity  on  the  open  boundary,  rf  and  u%  are  the 


/.  S/iu/man  etaL  /  Dynamics  of  Atmospheres  and  Oceans  61-62(2013)  35-45 


39 


MAY  JUN  JUL  AUG  SEP  MAY  JUN  JUL  AUO  SEP 


Fig.  3.  Sea  surface  height  (SSH)  anomalies  from  24  N  to  40°  N  as  computed  by  (left  panel)  the  regional  NCOM  CCS  model 
with  global  NCOM  boundary  conditions  and  (right  panel)  global  HYCOM.  The  observed  SSH  from  coastal  sea  level  observations 
at  Monterey  (36  36'  N)  and  San  Diego  (32°43'N)  are  Included  as  the  heavy  black  lines.  Time  series  at  the  same  locations 
sub-sampled  from  the  two  different  model  configurations  are  shown  as  the  red  (NCOM)  and  blue  (HYCOM)  lines  respectively. 

sea  surface  height  and  the  vertically  averaged  normal  component  of  velocity  from  larger  scale  model 
(NCOM  CCS  or  HYCOM).  The  open  boundary  condition  (1)  represents  a  radiation  condition  on  dif¬ 
ferences  between  the  Monterey  Bay  and  larger  scale  model  sea  level  elevations  and  transports.  The 
baroclinic  coupling  consists  of  using  the  vertical  structure  of  velocity,  temperature,  and  salinity  from 
the  larger  scale  model  in  specification  of  the  vertical  distributions  of  the  Monterey  Bay  model. 

3.  Experiments  design 

To  investigate  the  impact  of  remote  forcing,  model  resolution  and  bathymetry  on  coastal  model 
predictions  of  currents  on  the  shelf,  we  conducted  six  experiments  using  all  three  Monterey  model  set 
ups  (MBS1-3,  described  in  Section  2.2)  with  open  boundary  conditions  derived  from  two  larger-scale 
models  (NCOMCCS  and  HYCOM)  described  in  Section  23.  Attributes  of  six  runs  are  listed  in  Table  1. 


Tabic  1 

Monterey  Bay  model  runs. 


Runs 

OBCS 

Model 

Vertical  grid 

Bathymetry 

Horizontal  resolution 

1 

NCOM  CCS 

MBS1 

Sigma 

Smoothed 

13  km  x  23  km 

2 

NCOM  CCS 

MBS2 

Hybrid 

Non-smoothed 

13  km  x  23  km 

3 

NCOM  CCS 

MBS3 

Hybrid 

Non-smoot  hed 

03  km  x  0.7  km 

4 

HYCOM 

MBS1 

Sigma 

Smoothed 

13  km  x  23  km 

5 

HYCOM 

MBS2 

Hybrid 

Non-smoothed 

13  km  x  2.5  km 

6 

HYCOM 

MBS3 

Hybrid 

Non-smoothed 

03  km  x  0.7  km 

40 
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Run  1  was  conducted  with  the  MBS1  configuration  (Section  2.2).  The  run  was  initialized  on  27  July 
2006  00Z  from  the  NCOM  CCS  model,  and  run  until  1  September  2006  00Z  using  the  NCOM  CCS  as 
open  boundary  conditions. 

Run  2  was  conducted  with  the  MBS2  configuration  (Section  2.2).  The  run  also  was  initialized  on  27 
July  2006  00Z  from  the  NCOM  CCS  model,  and  run  until  1  September  2006  00Z  using  the  same  open 
boundary  conditions  as  in  Run  1 . 

Run  3  was  conducted  with  the  MBS3  configuration.  As  with  Runs  1  and  2.  Run  3  also  was  initialized 
on  27  July  2006  00Z  from  the  NCOM  CCS  model,  and  run  until  1  September  2006  00Z  using  NCOM  CCS 
as  open  boundary  conditions. 

Comparisons  of  Runs  1-3  demonstrate  the  impact  of  resolution  and  bathymetry  on  current  predic¬ 
tions  in  the  case  when  open  boundary  conditions  (NCOM  CCS)  underestimate  remote  forcing  in  the 
form  of  poleward  propagating  coastally-trapped  Kelvin  waves  (Section  23). 

Run  4-6  are  analogs  of  Runs  1-3  in  the  way  that  MBS1-MBS3  configurations  are  used  as  in  Runs 
1-3  (see  Table  1 ).  The  difference  is  that  Runs  4-6  were  initialized  from  the  global  HYCOM  model  on 
27  July  2006  00Z  and  were  run  with  the  open  boundary  conditions  from  HYCOM. 

Therefore,  comparisons  of  Runs  1-3  with  corresponding  Runs  4-6  demonstrate  the  impact  of 
remote  forcing  (in  the  form  of  poleward  propagating  coastally-trapped  Kelvin  waves)  on  the  model 
currents  predictions  on  the  shelf. 

The  magnitudes  of  complex  correlation  coefficients  and  angular  displacements  between  ADCP 
currents  and  model  currents  were  used  for  evaluating  the  above  model  runs. 

The  magnitude  R  of  the  complex  correlation  coefficient  and  the  angular  displacement  9  between 
the  ADCP  and  the  model  currents  fora  particular  depth  are  estimated  using  the  approach  outlined  in 
Kundu  (1976).  The  magnitude  R  is  estimated  as: 

R  =  V  Re2  +  Ini*  (2) 

where 

EtW  +  W 

\/£r(W)2  +  0'?)2)£rW2  +(^)2) 


_ Et(u°if  +  <W _ 

\/Er«u?)2+(l/?)2)Er((ur)2+(^)2) 


The  corresponding  angular  displacement  9,  which  is  also  called  the  phase  angle,  is  computed  according 
to: 


9  =  tan-1 


^((ufup  +  ifi^) 


(3) 


where  u^,  v {"  are  the  demeaned  east-west  and  north-south  model  velocity  components,  respectively, 
and  u°.  \Pt  are  the  demeaned  east-west  and  north-south  observed  velocity  components,  respectively. 
The  angular  displacement  9  gives  the  average  counterclockwise  angle  difference  between  model  and 
observed  velocity  vectors.  The  value  of  9  is  only  meaningful  if  R  is  significant.  Correlations  are  esti¬ 
mated  over  a  28-day  period  (August  2-30  of  2006)  using  hourly  model  and  observed  data  passing 
through  a  33-h  low-pass-filter.  Therefore,  the  actual  number  of  degrees  of  freedom  is  less  or  equal 
to  28  x  24/33  20.  With  20  degrees  of  freedom,  a  correlation  of  0.44  is  significant  at  95%  confidence 
level  (see,  for  example  Wilks.  1995).  Therefore,  correlations  less  than  0.44  should  be  considered  as 
insignificant. 


4.  Results 


Tables  2  and  3  present  complex  correlations  and  angular  displacements  between  observed  and 
model  currents.  Three  ADCP  bins  are  chosen:  near-surface  (12-1 6m),  mid-depth  (28-32,  48-52) 
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Table  2 

Complex  correlations  and  angular  displacements  between  model-predicted  and  observed  currents  at  ADCP1.  Statistically 
significant  correlations  and  corresponding  angular  displacements  are  In  bold. 


Depth 

12-16m 

28-32 

50-54 

R 

0 

R 

0 

R 

0 

Run  1 

0.1 

3  3 

0.12 

-135  3 

0.13 

123  3 

Run  2 

0.19 

-31 3 

0.13 

-66.0 

0.04 

1353 

Run  3 

0.11 

1093 

0.06 

-26.9 

0.14 

59.6 

Run  4 

039 

-19.8 

035 

-153 

036 

-7.6 

Run  5 

0.66 

-12.1 

0.79 

-6.7 

0.70 

-03 

Run  6 

037 

-15.6 

0.64 

-7  JO 

034 

-4.0 

Table  3 

Complex  correlations  and  angular  displacements  between  model-predicted  and  observed  currents  at  ADCP2.  Statistically 
significant  correlations  and  corresponding  angular  displacements  are  In  bold. 


Depth 

12-16m 

48-52 

88-92 

R 

0 

R 

0 

R 

0 

Run  1 

032 

175.5 

0.28 

1603 

0.08 

1503 

Run  2 

038 

433 

031 

23.1 

034 

-183 

Run  3 

0.21 

1013 

031 

1423 

034 

53 

Run  4 

0.63 

-133 

0.62 

-133 

0.45 

43 

Run  5 

0.83 

-43 

0.77 

-  63 

0.64 

-3.4 

Run  6 

0.40 

-213 

0.4 

-323 

030 

-28.8 

and  near-bottom  (50-54, 88-92)  for  ADCPs  (1,2),  respectively.  In  Tables  2  and  3,  complex  correlations 
which  are  larger  than  the  significant  level  (0.44,  Section  3)  are  in  bold  together  with  correspond¬ 
ing  angular  displacements.  All  three  model  runs  coupled  to  NCOM  CCS  (with  underestimated  remote 
forcing)  show  correlations  below  the  significance  level. 

Therefore,  as  long  as  the  Monterey  Bay  model  was  coupled  to  NCOM  CCS  (with  underestimated 
remote  forcing),  the  increase  in  horizontal  resolution  or  better  representation  of  bathymetry  did  not 
raise  the  model  versus  observed  currents  correlation  above  the  significance  level.  At  the  same  time, 
coupling  to  HYCOM  (with  more  accurate  representation  of  the  remote  forcing)  significantly  increased 
the  correlation  between  observed  and  model  currents.  For  Run  4,  with  coarse  resolution  and  smoothed 
bathymetry,  correlations  are  higher  than  the  significance  level  at  all  considered  depths  at  ADCP2,  and 
at  the  middle  depth  at  ADCPI.The  model  Run  5  (with  coarse  resolution  and  unsmoothed  bathymetry) 
performed  the  best:  for  all  depths  at  ADCPs  1  and  2,  correlations  are  higher  than  the  significance 
level  and  angular  displacements  are  not  larger  than  1 2°.  Results  in  Tables  2  and  3  indicate  that  the 
most  critical  element  in  reproducing  currents  on  the  shelf  is  an  accurate  representation  of  the  remote 
forcing  in  the  form  of  the  coastally-trapped  Kelvin  type  waves.  Comparisons  of  Runs  4  and  5  indicate 
that  the  second  most  important  factor  is  an  accurate  representation  of  bathymetry.  With  unsmoothed 
bathymetry  but  the  same  horizontal  resolution.  Run  5  has  a  better  correlation  with  observed  currents 
than  Run  4.  Model  Run  6  (with  finer  horizontal  resolution  than  in  Runs  4  and  5,  and  unsmoothed 
bathymetry  as  in  Run  5),  performed  worse  in  comparison  to  Run  5  at  ADCP1  and  Run  4  at  ADCP2. 

Visual  comparisons  of  observed  and  model-predicted  currents  are  shown  in  Figs.  4  and  5.  As  indi¬ 
cated  in  Ramp  et  al.  (201 1 ),  the  observed  currents  at  both  ADCPs  consisted  of  a  strong  poleward  flow 
during  3-8  August.  The  event  was  nearly  barotropic  at  both  locations.  This  poleward  flow  coincides 
with  the  strong  positive  SSH  anomalies  in  tide  gauges  data,  and  in  HYCOM  (Fig.  3).  Runs  1-3,  coupled 
to  the  NCOM  CCS  (that  underestimates  the  poleward  flow)  shows  weaker  and  more  oriented  inshore 
(positive  values  in  U-component  of  velocity)  poleward  flow  than  in  observations,  especially  at  the 
ADCP2  location.  Runs  4-6,  coupled  to  HYCOM,  captured  the  event  in  much  better  visual  agreement 
with  observations  (Figs.  4  and  5).  The  best  visual  match  in  reproducing  the  strength  and  direction  of 
this  poleward  flow  is  achieved  by  Run  5.  In  Run  6  (finer  horizontal  resolution  than  in  Runs  4  and  5),  the 
poleward  flow  is  weaker  and  somewhat  extended  up  to  August  12  at  ADCP2,  especially  below  20  m 
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Fig.  4.  Observed  (top  row)  and  model-predicted  U  (west-east,  eastward  flow  is  positive,  left  column)  and  V  (south-north, 
northward  flow  Is  positive,  right  column)  components  of  currents  at  ADCP1. 


depth  (Fig.  5).  This  is  probably  one  of  the  reasons  for  low  correlations  (below  the  significant  level) 
between  Run  6  and  ADCP2  currents  (Table  3). 

According  toRampetal.(2011)and  Figs.  4  and  5,  the  poleward  flow  was  followed  by  an  equatorward 
flow  (during  9-17  August)  at  mooring  locations.  Run  5  reproduced  this  reversal  best  in  timing  and 
intensity  (especially  at  ADCP2),  while  Run  6  shows  a  delay  of  the  reversal  below  20  m  depth  (almost  4 
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RUN  2 


RUN  3 


RUN  4 
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8  12  16  20  24  28 


8  12  16  20  24  28 


AUGUST  2006  (DAY)  vs.  DEPTH  (m) 


Fig.  5.  Observed  (top  row)  and  model-predicted  U  (west-east,  eastward  flow  is  positive,  left  column)  and  V  (south-north, 
northward  flow  is  positive,  right  column)  components  of  currents  at  ADCP2. 


days  delay).  All  Runs  1  -3  (coupled  to  NCOM  CCS)  missed  this  transition  from  poleward  to  equatorward 
flow,  especially  at  ADCP1. 

For  the  remainder  of  the  considered  record  (18-30  August),  observations  show  a  transition  from 
equatorward  flow  to  the  weak  poleward  flow  at  ADCP1  and  below  30  m  depth  at  ADCP2.  Model  Runs 
4-6  reproduced  this  weaker  poleward  flow  in  subsurface,  while  showing  equatorward  flow  at  surface 
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at  ADCP1  (up  to  10  m  depth),  and  deeper  (up  to  40  m  depth)  at  ADCP2.  Run  6  shows  a  delay  in  currents 
reversal  from  equatorward  to  poleward. 

5.  Conclusions  and  discussions 

Impacts  of  remote  forcing,  model  resolution  and  bathymetry  on  currents  predictions  at  two  ADCP 
sites  located  on  the  shelf  of  the  Monterey  Bay  area  are  investigated.  The  Monterey  Bay  area  model 
configurations  which  differ  in  horizontal  resolution,  vertical  grid  and  bathymetry  representation 
(smoothed  and  unsmoothed)  are  used.  Remote  forcing  (in  the  form  of  poleward  propagating  along 
the  coast  coastally-trapped  Kelvin  type  waves)  introduced  through  the  coupling  of  the  Monterey  Bay 
model  configurations  to  the  larger  scale  models:  the  global  HYCOM  with  accurate  representation  of 
remote  forcing,  and  the  regional  California  Current  model  (NCOM  CCS)  with  underestimated  remote 
forcing. 

Comparisons  of  complex  correlations  and  angular  displacements  between  observed  and  model 
currents  show  that  all  model  runs  coupled  to  NCOM  CCS  (with  underestimated  remote  forcing)  show 
correlations  below  the  significance  level.  At  the  same  time,  coupling  to  HYCOM  (with  accurate  repre¬ 
sentation  of  the  remote  forcing)  significantly  increased  the  correlation  between  observed  and  model 
currents.  Even  the  model  run  (Run  4  in  Table  1)  with  coarse  resolution  and  smoothed  bathymetry 
showed  higher  than  the  significant  level  correlations  at  all  considered  depths  at  ADCP2,  and  at  the 
middle  depth  at  ADCP1.  These  results  indicate  that  the  most  critical  element  in  reproducing  currents 
on  the  shelf  is  an  accurate  representation  of  the  remote  forcing  in  the  form  of  the  coastally-trapped 
Kelvin  type  waves.  Also,  comparisons  of  complex  correlations  have  shown  that  an  accurate  represen¬ 
tation  of  bathymetry  is  the  second  most  critical  factor  to  remote  forcing  in  reproducing  currents  on 
the  shelf. 

Visual  comparisons  of  observed  and  model-predicted  currents  support  the  above  conclusions.  Runs 
coupled  with  the  HYCOM  reproduced  much  better  the  intensity  and  direction  of  the  observed  strong 
poleward  flow  (3-8  August,  Figs.  4  and  5).  This  poleward  flow  coincides  with  the  strong  positive  SSH 
anomal  ies  presented  in  tide  gauges  observations  and  the  HYCOM  (Fig.  3).  Run  5,  coupled  to  the  HYCOM 
and  with  unsmoothed  bathymetry,  reproduced  the  best  timing  of  the  observed  transition  from  the 
poleward  flow  to  equatorward  flow  (observed  during  9-17  August),  and  the  transition  back  to  the 
weak  poleward  flow  for  the  reminder  of  the  considered  record  (18-30  August). 

As  shown  in  Section  4,  the  impact  of  the  horizontal  resolution  on  the  model  predictions  left  mixed 
results.  The  model  run  with  the  finer  horizontal  resolution  and  unsmoothed  bathymetry  (Run  6. 
MBS3  model  configuration)  performed  worse  in  comparison  to  the  run  with  coarse  resolution  and 
unsmoothed  bathymetry  (Run  5.  MBS2  model  configuration)  at  both  ADCP  locations.  Also,  Run  6 
demonstrated  worse  predictions  at  the  ADCP2  location  in  comparison  to  Run  4  with  coarser  resolu¬ 
tion  and  smoothed  bathymetry  (MBS1  model  configuration).  The  decrease  in  predictive  skill  as  we 
approach  a  very  fine  grid  resolution  might  be  the  result  of  such  factors  as  the  deficiencies  in  physi¬ 
cally  based  representation  of  dissipation  rates,  drag  laws,  sub-grid-scale  parameterizations,  truncation 
errors,  as  well  as  by  approaching  the  limit  of  hydrostatic  assumption.  Another  possible  explanation 
might  be  the  ratio  of  horizontal  grid  resolutions  between  the  Monterey  Bay  model  and  the  larger 
scale  model.  The  HYCOM  model  resolution  is  around  9  km,  while  the  MBS3  configuration  has  around 
1  -1.5  km  resolution  on  open  boundaries. This  provides  approximately  1  to  9  or  1  to  6  ratios  in  one  way 
coupling  on  open  boundaries.  These  high  ratios  of  grid  resolutions  between  larger  scale  versus  local, 
finer  resolution  models  in  one  way  coupling  probably  introduce  artificial  time  lags  in  propagating 
information  from  the  larger  scale  model  to  the  finer  resolution  Monterey  Bay  model  (MBS3  configu¬ 
ration).  This  might  introduce  time  lags  in  reproducing  reversals  from  poleward  to  equatorward  flows 
and  back,  which  are  discussed  in  Section  4.  The  above  issues  with  the  MBS3  configuration  are  planned 
to  be  investigated  in  our  future  research. 
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